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The matching conditions at the interface between two non-mixed fluids at rest are obtained 
directly using the equation of movement of the whole media. This is a non-usual point of view 
in hydrodynamics courses and our aim is to fix ideas about the intrinsic information contained 
in the matching conditions, on fluids in this case. Afterward, it is analyzed the polarization of 
the normal modes at the interface and it is shown that this information can be achieved through 
a physical analysis and reinforced later by the matching conditions. A detailed analysis of the 
matching conditions is given to understand the role that plays the continuity of the stress tensor 
through the interface on the physics of the surface particle movement. The main importance of the 
viscosity of each medium is deduced. 



En el presente trabajo se utilizan las ecuaciones de movimiento de todo el sistema compuesto 
por dos fluidos no miscibles en reposo relativo para obtener las condiciones de empalme en la 
intercara. Este procedimiento es inusual y el objetivo de hacerlo en esa forma es fijar ideas sobre la 
informacion que esta contenida en dichas condiciones de empalme, en este caso en fluidos. Posterior- 
mente es analizada la polarizacion de los modos normales en la interface y se demuestra que estas 
caracten'sticas pueden ser obtenidas directamente de analisis fi'sico de las condiciones de empalme. 
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I. INTRODUCTION 



The boundary conditions constitute the key feature in 
any theory of surface waves. It is through them that 
one introduces in the analysis the physical consequences 
of specific surface effects, such as local changes of mass 
or stresses, thus going beyond the simple case in which 
one merely matches two semi-infinite bulk media at an 
interface (which, in particular, can be a free surface). 

In the last years, increased attention has been paid 
to the properties of capillary waves by physicists and 
chemists Jl-§]. Ripples represent one of the few cases 
in which the relation between the dynamical properties 
of a surface and liquid flows can be predicted completely. 
The study of capillary ripples has clarified which proper- 
ties of a liquid surface determine the surface's resistance 
against deformation. 

The boundary conditions for the stress at the interface 
are derived from the principle that the forces acting upon 
such an "interfacial" element result not only from viscous 
stresses in the liquid but also from stresses existing in the 



deformed interface. The cause of the difference is that 
an interface, unlike a three dimensional liquid, can not 
enjoy the property of incompressibility Q. Work done 
on an element of liquid is partially degrade into heat by 
viscous friction and partially transformed into kinetic en- 
ergy which is transmitted to adjoining elements. Work 
done on an interfacial clement leads, at least partially, to 
an increment of the surface potential energy. It is this 
potential energy of the deformed interface which enables 
the whole system, including the interface, to carry out 
an oscillatory motion. 

In this article, first, it will be shown a non-traditional 
way to obtain the matching conditions at the interface be- 
tween two non-mixed fluid at rest, considering the equa- 
tion of motion of the whole media directly. Usually, uni- 
versity courses do not use this approach and state the 
boundary conditions from outside the constitutive equa- 
tions governing the studied problem. It is important from 
a pedagogical point of view to evidence that in fact, the 
matching conditions at the interfaces arc contained, al- 
most in all cases, in the equation of movement for the 
whole system taken as the composition of all media. The 
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other cases arise when the fluid interfaces have intrinsic 
properties not included on the equations of the media in 
that cases mentioned before. 

This point of view was already used in [p|-^3| to elab- 
orate a formalism, based on Surface Green Functions, 
which establish an isomorphism between solids and flu- 
ids related to the interface normal modes. The aim of this 
paper is to emphasize the fact that this way to establish 
the matching conditions allows us to introduce the phys- 
ical characteristics related with the boundary conditions 
in a firmer floating. 



II. EQUATIONS FOR THE BULK IN A VISCOUS 
INCOMPRESSIBLE FLUID. 

To achieve a system of equations which describes the 
oscillatory motion of a viscous incompressible fluid, the 
starting point is the linearized Navier-Stokes equation for 
a viscous fluid ifESl, which reads: 
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where p is the density of equilibrium, Vi the components 
of the fluid particle, F t ex ^ are the external forces and ry 
is the stress tensor which, for a viscous incompressible 
fluid, has the following form 
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Here p are the small variations of pressure, rj is the 
viscosity parameter and Sij is the unit matrix elements. 
The subscripts i and j take the values y and z identically. 
In (Q) the sum over repeated subscripts is understood. It 
had been assumed that the system is symmetric with re- 
spect to the x direction. 

Putting (0) in (fy) it is obtained 
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where we have neglected the external forces. 

Also the continuity equation is needed, which expresses 
that the volume of an element of the incompressible fluid 
does not change during the motion. It has the form (13 
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The first step is to obtain the equation of motion and 
its solution of the each medium taken as infinite. So, 



the parameters as density and viscosity are considered 
constants in the whole medium and it leads to transform 
equations (§) and (@) to: 
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where V 2 = (d 2 /dy 2 + d 2 /dz 2 ). 

The solution of the system (g), and (0) written as a 
vector field velocity, can be putted as the sum of an irro- 
tational field (related with the longitudinal mode) and a 
divergence free field (related with the transverse modes) 
|, i.e.: 



which satisfy: 



V = V 1 + V 2 
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Any irrotational field is characterized by a scalar func- 
tion, the "potential" function ip(y, z, t), such that 
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and the divergence free field can be described by a vector 
function Jl4| , the 11 stream" or vorticity function ip(y, z, t), 
such that 
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The velocity components can thus be written in terms 
of the potential and the stream functions, as: 



1 = Jlf_ _ 

y dy dz 

Vz = _d 1 + d± 
dz dy 



(13) 
(14) 



Substitution of eqns. ([13]) and ( pT| ) in the continuity 
condition (||) gives 

VV = (15) 
and substitution into eqns. @ and (^) leads to 
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Equations (|16|) and (17) are simultaneously satisfied if 
one considers: 
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(18) 
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The constants C\ and C 2 are obtained from the condi- 
tion at zero flow, and it gives raise to C\ = p and C 2 — 
where p a is the reference atmospheric pressure. 

The solution of equations (|l5|), ( fl8| ) and (|l^) are look- 
ing for as the following 



< P (v,z,t) = $(z)e i ( K v-ut) 



(20) 
(21) 



where the parameters k and u> are the wavevector and 
frequency of the wave respectively. 

Substituting eqns. @ and @ in @ and © gives 
the z-dependence of the functions p and ip, which satisfy: 
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with q 2 = k 2 — ipui/r] 

Equations ( p2| ) and ( |23| ) lead to the solutions: 
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and in combination with eqns. ( |20| ) and (|2l]) they give a 
solution of the form: 



<p( y , z, t) = Id e KZ +c 2 e- KZ ) e*( K v - ut ) (26) 
1%, z, t) = (c 3 +c 4 e~%«) e*(^ - ut ) (27) 



where Ci, C2, C3 and C4 are constants to be determined 
by boundary and matching conditions. 

Then, the expressions for the varying velocity com- 
ponents and the pressure are finally obtained by substi- 
tution of eqns. @ and @ into @, (Q) and @, 
respectively. 



III. INTERFACE PROBLEM: MATCHING 
CONDITIONS. 

Now we will match the two media. We consider a sur- 
face which at rest coincides with the plane z = and 
it separates medium Mi at z < from medium M 2 at 
z > 0. Each one is viscous and incompressible. 

First of all, any solution has to fulfill the continuity of 
the velocity field across the surface according to 
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Superscript 1 denotes medium Mi and superscript 2 
denotes medium M 2 . The rest of the matching conditions 
at the interface are derived from the system of equations 
which govern the whole system. This is no usually done 
in the normal program courses at the Universities and we 
consider that it is important to state that the matching 
conditions are, almost in all of the cases, contained in 
the equation of movement for the whole system taken as 
the composition of each one. These are eqns. (||) and (Q) 
where the parameters r\ and p are constant, but taking 
different values on each medium. Integrating these equa- 
tions through the surface about z = from — e to +e and 
later taking e — > 0, it is obtained from eq. (^) 
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and from eq. 
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where p 1 is the jump due to the surface tension according 
with the Laplace Law |l4[ . The other terms in eqns. (||) 
and (0) vanish when e —> because they are continuous 
or have a finite jump at z = 0. 

It is easily seen according to eq. (||) that eqns. ( |30| ) 
and ([H]) are the conditions of the continuity of the stress 
tensor components through the interface, as expected. 

On the other hand, the boundary conditions of the 
problem are regularity at z — > ±00, which leads to elim- 
inate 4 of the 8 constants appearing in ( p6| ) and ( |27| ) for 
the two media. Using the resulting functions ip and ip for 
each medium in conditions (|2^)-([5l]) the following system 
is met: 
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to determine the remaining constants, with the defini- 
tions da = —rjiQti—iY, c?43 = 2i«77iq^ 1 — T, T = /ui, 
Qn = 2k 2 - ipiw/171 and Q t2 = 2k 2 - ip 2 uj/ri 2 . 

In deducing the expression in system (^2|) which comes 
from (HI) it was considered that 
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on the surface, but as we are dealing with the velocity 
field, then it is obtained that: 
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which finally leads to 
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according to the substitution d/dt — > — iui and d/dy — > 

The vanishing of the determinant of ( |32| ) leads to the 
equation for the dispersion relation of the existing surface 
modes. It can be written as: 

[(pi + p 2 )(piqt 2 + mti) - k (pi - P2) 2 ] + 

+ 7K 3 [pi(K-q t2 )+p 2 (K-q t i)] + 
+ 4k 3 (?72 - ??i) 2 (k - qti)(« - qt2) + 
+ AiK 2 uj(r]2 - m)(pi K ~ P2K - Piqt2 + Paqti) = (36) 

This is the dispersion relation for capillary waves for 
the interface of two viscous incompressible fluids. 

In order to obtain the constants Ci, C2, C3 and C4, an 
initial stimulus is needed according to an initial value 
problem Jlq ] but this method, simple at the beginning, 
becomes rather complicated later and it is not good for 
a quite general study. 

As was said on the introduction, the aim of this paper 
is not to get inside the dispersion relation of the nor- 
mal modes at the interface of two viscous fluids at rest. 
For a better study of this subject we recommend paper 
p6| . We get here to show a way of solution also different 
from the usually taken as only the velocity vector as a 
function of an " stream function" . From now on, we will 
put our attention on the matching conditions and we will 
show that more that a mathematical information of the 
matching can be found on it, but also the physics of the 
polarization can be deduced and how the interface moves 
in its oscillation. 

From eqns. © and © it can be obtained more 
information about the polarization of the modes on the 
interface. This will be done in the next section. 



mean that there are no longitudinal modes of oscillation 
on the interface when the fluids involved are incompress- 
ible. It had been shown in [|) that the interface, when 
oscillating, must be considered as a compressible one, be- 
cause precisely its change in area is responsable for the 
increasing of its potential energy and therefore, for its 
oscillation. It is important to state that this is a fun- 
damental argument to understand the movement of any 
interface in hydrodynamics. 

Nevertheless, now it can be shown that in spite of the 
compressibility of the interface, there does not exist pure 
longitudinal mode if we are dealing with incompressible 
media. Let us demonstrate this. 

If a point y on the interface is considered moving with 
velocity, say Vy in the y-axis direction, then, according 
to the continuity of velocity, the point (y ,— e) in Mi and 
the point (y ,+e) in Mi must have the same velocity Vy 
if e is small enough. As the interface is compressible, 
at the point y\ near enough y a the velocity can be, for 
instance, Vy-y different in general from Vy but as the me- 
dia are incompressible, at the point (j/i,— e) and (yi,+e) 
the velocity must be Vy . See Fig. 1. This is not in 
conformity with the continuity of the velocity through 
the surface and hence the pure longitudinal mode is not 
possible and only the TN mode seems to be valid when 
the media are incompressible. 

After these considerations during the above demon- 
stration, the student can keep the idea that the interface 
oscillations can only occur in the z-axis. This is the ac- 
curate moment to show to the student that things not 
always are as they apparently seem to be, because that 
assumption does not take into account the different prop- 
erties of each medium, whose response depends on its 
fundamental parameters, as density and viscosity, which 
are different for each medium. Hence, it is evident that it 
must be analyzed, precisely, the interface matching con- 
ditions. 

It is useful to compare and to support the previous 
qualitative analysis with a quantitative and more pro- 
found one regarding the interface matching conditions. 

Recalling carefully eqns. (|0|) and (|3l]) and supposing 
that such a wave propagates in y direction with move- 
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ment only in z direction (TN mode), then Vy — v y 
and eq. (p5C|) becomes 
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IV. MATCHING CONDITIONS AND 
POLARIZATION. 
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There are two possible modes on the fluid: one in which 
the fluid particle moves in the direction of the wave prop- 
agation called longitudinal with notation L(Vy,0) and an- 
other transverse to the direction of the wave propagation 
and normal to the interface with notation TN(0,Vz). 

The longitudinal mode is related to the fluid compress- 
ibility because this motion of the fluid particles is only 
possible when its volume changes |l4|] . This analysis also 
holds when the mode is on the interface, but this docs not 



It is known that V z is continuous along the interface for 
all points. Then, the derivative with respect to y is also 
the same in both hands of (|37j ) and this expression only 
holds if 771 = 772, i.e., if the media have the same viscosity. 
It does not mean for the interface to disappear because 
the density of each medium can be different. Then, if 
the viscosities of the media have not the same value, the 
velocity component Vy along the interface must be differ- 
ent from zero to compensate the inequality (p7[) and to 
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fulfill the continuity of the stress tensor in the y direction 
given by eq. (|30| ) yielding to a component of movement 
along y direction. This mode will be called Sagittal mode 
orS(V y ,V z ). 

The above analysis was done for the general case. Now 
we are able to take the particular case in which one of the 
media is vacuum, for instance, M%, with 772 = 0. Then, 
condition (Esol) becomes 



'6V§> 
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and it can be seen that also Vy must be non zero on the 
surface to hold eq. ( |3^ ) with the corresponding Sagittal 
polarization movement. 

With this analysis on the conditions of stress compo- 
nent continuity in y direction along the interface, it can 
be seen that if the two media are viscous (at least one of 
them), the fluid particle of the interface moves in a Sagit- 
tal mode which combines movement in both directions: 
along the wave propagation in y direction, and normal to 
the interface in z direction. 

Then, it is qualitatively clear that the viscosity of each 
media plays a fundamental role in the coupling of modes 
even for incompressible fluids. In spite of that, it could be 
a mistake to say that the modes decouple if the viscosities 
are equal. It should not be forgotten that eq. (|3l]) is also 
important in the characterization of the interface particle 
behaviour and it includes the pressure on each side of the 
surface. According to eq. (jig), the pressure is associated 
with the longitudinal mode and the inertial effect of the 
fluid particle according with the density of the media. 
Then, it contains the information of each components of 
the velocity and also of the density and according to eq. 
( p6| ) the pressure has a jump through the interface. This 
result, in combination with the analysis of eq. ( |30| ) make 
difficult to understand the role played by the densities of 
each media on the surface polarization movement, and it 
can not be reached from this only analysis. This point is 
still a matter of investigation. 



V. CONCLUSIONS 

The present work is an attempt to give an example, us- 
ing the hydrodynamics, of how the study of the interface 
matching conditions allows us to make a plentiful and 
rich in details discussion. Moreover, of how the inter- 
face matching conditions content a sufficient information 
to conclude that the interface oscillation must be with 



a Sagittal mode and not with neither a pure longitudi- 
nal, nor a pure transversal one. This movement has been 
shown to be close related to the physical properties of 
the media such as viscosities and that fact allows us to 
establish rigorously that those are the parameter which 
characterize the interface movement and the response of 
each medium to an stimulus coming from the other one. 

It was seen that viscosity is the main parameter in 
the coupling of the two modes to achieve a Sagittal one, 
nevertheless within the framework of this formalism it is 
difficult to determine the role of viscosity and of the den- 
sity ratios in the coupling of modes. This aspect needs 
further investigation. 
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FIG. 1. Relation of velocities on the interface between two 
viscous fluids. 
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